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ABSTRACT 

o 

• Context. The outflow driven by the low-mass class protostar LI 157 is the prototype of the so-called chemically active outflows. The bright 
bowshock B 1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and 
chemical processes are strongly interdependent. 
JL Aims. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key prograin 
>| "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the 

excitation conditions in the shock region. 
^ Methods. The CO 5-4 and 0-H2O lio - loi lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 
< • lb spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window 

help estabUsh the origin of the emission detected in the main-beam of HUT, and the physical conditions in the shock. 
^Sj Results. Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission 
^ arises from two regions with distinct physical conditions : an extended, warm (100 K), dense (3 x 10' cm"') component at low-velocity, which 
dominates the water line flux in Band 1 ; a secondary component in a small region of B 1 (a few arcsec) associated with high- velocity, hot (> 400 K) 
gas of moderate density ((1.0 - 3.0) x lO'* cm"'), which appears to dominate the flux of the water line at 179 /jm observed with PACS. The water 
fT^ abundance is enhanced by two orders of magnitude between the low- and the high- velocity component, from 8 x 10"' up to 8 x 10"' . The properties 
•T-H of the high-velocity component agree well with the predictions of steady-state C-shock models. 

^ Key words. ISM: individual objects: LI 157 — ISM: molecules — stars: formation 

o 
o 

". 1 1 . Introduction several orders of magnitude, as reported for various molecular 

^ species in "chemically active" outflows (Bachiller et al. 2001). 

k> Shocks m protostellar outflows play a crucial role m the molec- Conversely, the magnetic field and the ionization fraction play 

^ ular cloud evolution and star formation by transferrmg momen- ^„ important role in controlling the size and the temperature of 

H turn and energy back to flie ambient medium. There is mount- t^e ion-neutral drift zone. Because of the interplay between the 

ing evidence that these shocks often involve a magnetic precur- dynamics and chemistry the physics of MHD shocks requires a 

sor where ionic and neutral species are kinematically decoupled, comprehensive picture of bofli the gas and dust physical condi- 

Magneto-hydrodynamical (MHD) shocks are important not only ^^^^^ ^jj^ compressed region itself, 
for the cloud dynamics, but also for the chemical evolution 

through temperature and density changes, which favors the acti- Along with H2, H2O and CO are two key-molecules pre- 

vation of endothermic reactions, ionization, and dust destruction dieted to dominate the cooling of MHD shocks (Kaufman & 

through sputtering and shattering in the ion neutral drift zone. Neufeld 1996). The abundance of H2O in protosteUar regions 

These various processes lead to abundance enhancements up to can be greatly enhanced in shocks, even of moderate velocity. 

This occurs both from the sputtering of frozen water from grain 

* Herschel is an ESA space observatory with science instruments mantles and tiirough high-temperature sensitive reactions in the 

provided by European-led Principal Investigator consortia and with im- gas phase (EHtzur & de Jong, 1978; EHtzur & Watson, 1978; 

portant participation from NASA. Bergin et al. 1998). Multi-transition observations of these two 
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molecules therefore serve as good probes of shock regions with 
various excitation conditions, and can be used to set stringent 
constraints on MHD shock models (Flower & Pineau des Forets, 
2010). 

The heterodyne instrument HIFI onboard Herschel allows us 
to study with unprecedented sensitivity the chemical and dynam- 
ical evolution of protostellar shocks, at spectral and angular res- 
olutions comparable to the best ground-based single-dish tele- 
scopes. This is the main goal of the spectral survey of LI 157-B 1, 
carried out in the guaranteed time key-project CHESS. 

The source L1157-mm is a low-mass Class protostar lo- 
cated at a distance estimated between 250 pc (Looney et al. 
1997) and 440 pc (Viotti 1969). It drives a spectacular bipolar 
outflow, which has been studied in detail at millimeter and far- 
infrared wavelengths. Mapping of the southern lobe of LI 157 
with the Plateau de Bure Interferometer (PdBI) reveals two limb- 
brightened cavities (Gueth et al. 1996), each one terminated by a 
strong bow shock, dubbed "Bl" and "B2" respectively (Fig. 1), 
which are likely the result of episodic ejection events in a pre- 
cessing, highly collimated jet. The spatial and kinematical struc- 
ture of Bl has been modelled in great detail by various au- 
thors, making it the archetype of protostellar bowshocks in low- 
mass star-forming regions and the testbed of MHD shock models 
(Gusdorfetal. 2008a,b). 

Here, we report on the emission lines of CO and H2O de- 
tected in the low-frequency band of HIFI in the course of the 
CHESS spectral survey. From comparison with complementary 
observations, we discuss the origin of the emission, and based on 
a simple modelling of the source, we derive the water abundance 
in the shock region. 

2. Observations and results 

A fun coverage of the band lb at 20'' 39"' 10.2' +68°01'10.5" 
(J2000) in the bowshock B 1 was carried out with the HIFI het- 
erodyne instrument (de Graauw et al. 2010) on board of the 
Herschel Space Observatory (Pilbratt et al. 2010) during the 
Performance Verification phase on 2009 August 1. The corre- 
sponding dataset is OBS_1342181 160. The HIFI band lb (from 
555.4 to 636.2 GHz) was covered in double beam switching. 
Both polarizations (H and V) were observed simultaneously. The 
receiver was tuned in double sideband (DSB) with a total inte- 
gration time of 140 minutes. In order to obtain the best possible 
data reconstruction, the survey was acquired with a degree of re- 
dundancy of 4. The Wide Band Spectrometer (WBS) was used 
as spectrometer, providing a frequency resolution of 0.5 MHz. 

The data were processed with the ESA-supported package 
HIPE (Ott et al. 2009). Fits file from level 2 were then cre- 
ated and transformed into GILDAS formaf] for baseline sub- 
traction and subsequent sideband deconvolution. The spectral 
resolution was degraded to IMHz in the final single sideband 
(SSB) dataset. The calibration for each receiver (H and V) is bet- 
ter than 2-3%. The relative calibration between both receivers is 
also rather good, with a difference in intensity of about 4%. The 
overall calibration uncertainty is about 7%, except for the strong 
CO line present in the band (see below). 

Two strong lines dominate over the molecular transitions de- 
tected in the spectral band: the fundamental line of water in its 
ortho state 0-H2O lio - loi at 556.936069 GHz and the CO 5- 
4 fine at 576.276905 GHz (Fig. 1). The final rms noise is 13 
mK. We adopted the theoretical telescope main-beam efficiency 
Tjinb - 0.723, and a main-beam size of 39" (HPFW) in the whole 
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Fig.l. (left) Southern outflow lobe of LI 157 in CO 1-0 
(greyscale and black contours) and in SiO 2-1 (white contours) 
as observed at the PdBI (Gueth et al. 1996,98). A black square 
marks the nominal position of bowshock L1157-B1 observed 
with HIFI. The HIFI main-beam is represented with a black cir- 
cle, (right) Panel of the CO 5-4 (bottom) and 0-H2O lio - loi 
(centre) line spectra obtained in band lb of HIFI. For both lines, 
we show (dashed) a magnified and spectrally smoothed view of 
the emission. Intensities are expressed in units of antenna tem- 
perature, (fo/9) H2O spectrum with fitted low-velocity component 
(LVC), high-velocity component (HVC), absorption feature and 
summed fitted spectrum. 



band. Unless indicated, intensities are expressed in units of an- 
tenna temperature Ta. 

The CO 5-4 transition is detected with an intensity of 9 K 
(Ta) and a linewidth of 5 km s"'. We notice a weak absorption 
feature in the line profile in the redshifted gas over a wide ve- 
locity range, which may partly arise from cloud contamination 
in the reference position. The intensity in the blue wing of the 
CO line differs by as much as 20% between both polarizations. 
This effect is not observed towards the H2O line. Its origin is not 
understood at the moment. 

The fundamental 0-H2O lio - loi line is detected with an 
intensity of 0.9 K (Ta) at the peak. It is characterized by a 
broad linewidth ^ 15kms The line displays an absorption 
dip at visi = -1-2.9 km s"' and a broad redshifted wing extend- 
ing up to H-8kms The broad linewidth of the H2O spec- 
trum could be fit with three Gaussian velocity components, a 
low-velocity, a high-velocity, and an absorption component. The 
low (high) velocity component peaks at v^r = -0.58 km s"' 
(-7.86 km s '); the linewidth and peak intensity derived from 
the fit are 9.56 km s-i and 0.77 K (13.72 km s"' and 0.29 K), re- 
spectively. The absorption component was fit by a narrow line 
Gaussian (AV = 1.38 kms"') of amplitude -0.48 K centered at 
visr - +2.9 km s ' . The fit of the individual components and the 
resulting fit to the water spectrum is displayed in Fig. 1 . 

Overall, the H2O and CO emission are detected in the same 
velocity range. The high sensitivity of the HIFI observations 
permits the detection of emission from the entrained gas up to 
VhT = -30kms"', i.e. about lOkms"' higher than was pre- 
viously known from ground based observations. However, line 
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Fig. 2. {left) Velocity -integrated CO 6-5 emission maps of the 
low- and high-velocity components. LVC (HVC) emission is 
represented in greyscale and thin dashed contours (white con- 
tours); first contour and contour interval are 3cr and Icr (10% of 
the peak flux), respectively, Contours at half-power are drawn in 
thick, (right) Same for SiO 2-1 observed at the Pdbl. The fflFT 
main-beam is represented with a black circle. 



profiles differ noticeably and the ratio of the H2O/CO5 - 4 line 
intensities increases with increasing velocities from about 0.2 in 
the ambient gas up to 0.9 at visr = -25 km s'' (Fig. 3). 

All the other molecular tracers detected in the HIFI band 
show a pronounced break in the line profile at visr ~ -7.2 km s"' 
(Codella et al. this volume). This is also observed in the CO 
6-5 spectrum of Bl obtained by us at the CSO, as part of com- 
plementary observations to help analyse the CHESS data. The 
maps of the whole southern lobe of LI 157 in CO 3-2 and 6-5 
obtained at the CSO in June 2009, with 24" and 14.5" respec- 
tively, will be discussed in detail in a forthcoming paper (Lefloch 
et al. 2010, in prep.). 

Below, we define the region with visr < -7.25 km s~' as the 
high-velocity component, hereafter HVC (see also Codella et 
al), and the region with visr > -7.25 km s"') as the low-velocity 
component (LVC). As we discuss below, these two velocity com- 
ponents are characterized by different spatial extents and excita- 
tion conditions. 



3. Discussion 

3. 7. Ongrn of the emission 

Due to its relatively high energy above the ground state (£■„,, = 
1 16 K) the CO 6-5 transition is a good probe of the warm regions 
where H2O can evaporate from grain mantles and be released in 
the gas phase. SiO 2-1, observed at the PdBl at 2.5" resolution 
is a particularly good tracer of shocks strong enough to release 
refractory elements in the gas phase, because it is usually unde- 
tected in the cold, quiescent molecular gas. 

The overall SiO emission is strongly peaked at the position 
of Bl, which appears as a region of ^ 15" size located at the 
apex of the cavity. Interferometric maps of the southern lobe 
(Figs. 1-2) reveal extended emission along the eastem wall of 
the cavity (the low- velocity wing of the bow) and downstream of 
Bl, at velocities close to systemic, both blue and red. By com- 
paring these data with IRAM 30m observations (Bachiller et al. 
2001), we checked that unhke the HVC, a fraction of the flux 
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Fig. 3. (bottom left) Comparison of the 0-H2O lio- loi line 
profile with the SiO 2-1 emission observed at the PdBl, 
averaged over the HlFl beam, (top left) Variations of the 
SiO 2-1 / H2O line ratio as a function of velocity, (bottom 
right) Comparison of the 0-H2O lio - loi and CO 5-4 line pro- 
files, (top right) Variations of the H2O /CO 5-4 line ratio as a 
function of velocity, smoothed to a resolution of 2 km s"^ . 



emitted in the LVC is actually missed in the PdBl data, which 
is direct evidence for extended emission. This is consistent with 
the CO 6-5 data (Fig. 2). The low- velocity gas emission is lo- 
cated in the wake of Bl, reaching ~ 40" North from the apex. 
The area of the LVC amounts to « 1/3 of the HIFI beam (see 
Fig. 2). Interestingly, the PACS map of the H2O 179 yum line 
reveals large-scale emission, spatially coinciding with SiO 2-1 
emission in the outflow (Nisini et al. 2010). 

In any case, there is definitely much less molecular gas emis- 
sion associated with the western wall of the cavity (Benedettini 
et al. 2007). We therefore expect an asymmetry in the H2O spa- 
tial distribution, similar to that observed in many other tracers 
such as CS or HCN, as shown by the PACS map of the 179 jum 
H2O line (Nisini et al. 2010). 

We note an excellent agreement between the H2O and the 
low-excitation SiO 2-1 line profiles (Fig. 3) in the high- velocity 
range, with a constant SiO 2-1 / H2O line ratio » 0.8 between -20 
and -7kms This has important implications. First, this con- 
stant ratio in the range of the HVC suggests that both emissions 
most likely arise from the same region and that the emissions are 
optically thick. In that case, the low intensities measured in the 
high- velocity component (a few tenths of K) point to a small size 
extent. This is direct evidence that the H2O emission detected 
fiUs only partly the HIFI beam. Indeed, the bulk of the SiO HVC 
originates from a small region of 4" X 12" in B 1 (Fig. 2), cor- 
responding to a filling factor fF ~ 0.03 in the HIFI main-beam. 
Second, if silicon comes from grain erosion, the SiO profile is 
predicted to be much narrower than H2O because it takes a long 
time for Si to oxidize into SiO, so SiO comes only from the cold 
postshock, as discussed by Gusdorf et al. (2008b). The similarity 
of the SiO and H2O line profiles suggests that SiO forms more 
extensively in the shock than predicted by oxidation of sputtered 
Si atoms. As it can be released in the gas phase even at low ve- 
locities in the shock, SiO is present in the gas phase over the full 
width of the shock wave. 

In summary, we find strong observational evidence that the 
emission from the HVC and LVC arises from regions of differ- 
ent physical extent. The size of the HVC appears definitely much 
lower than the LVC (// ^ 0.03 and 0.3, respectively). It is true 
however that the present determinations are uncertain. HIFI ob- 
servations of the high-excitation lines of CO and H2O will make 
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Table 1. Observed and fitted parameters of the H2O and CO 5-4 lines and prediction for the 179 /im H2O line flux. The CO fluxes 
in the HVC and LVC are integrated in the velocity intervals [-30;-7.25] and [-7.25; +1 1.0], respectively. The H2O fluxes are derived 
from a multiple Gaussian fit to the line profile. Intensities are expressed in units of antenna temperature (Ta) K km s"' . 





H2O 557GHz 


CO 5-4 


CO 6-5' 


CO 3-2' 


Size 


N(CO)" 




ji 


X(H20)' 


F(179yum) 




(Kkms-i) 


(Kkms-i) 


(Kkms-i) 


Kkms-' 


(") 


( cm"^) 


( cm-3) 


(K) 


( cm-2) 


(W cm-2) 


LVC 


7.83 


45.4 


29.2 


40.4 


25 


8.0(16) 


(1.0-3.0)(5) 


100 


0.8(-6) 


4.2(-20) 


HVC 


4.28 


3.98 


1.34 


3.11 


7 


5.0(16) 


(1.0-3.0)(4) 


400 


0.8(-4) 


7.1 (-20) 



' From CSO observations smoothed down to the resolution of the HIFI observations. 

^ Determined from LVG analysis of the CO emission. 

^ From comparison of LVG-derived N(o - H2O) with N(CO), assuming a water OPR of 3 and an abundance [CO]/[H2] = 10 
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Fig. 4. Best-fit solution to the LVG modelling of CO line tem- 
peratures in the 3-2, 5-4 and 6-5 transitions for both HVC (left) 
and LVC (right) components, respectively. The contour of the 
observed CO 5-4 integrated line area (TAv), corrected for main- 
beam dilution, is drawn in a dashed line, as well as the 5-4/3-2 
and 6-5/3-2 line ratios (solid and dotted lines, respectively). Thin 
lines deUneate the uncertainties in the observed values. 



it possible to better estabhsh appropriate filling factors for these 
components. 

3.2. Physical conditions 

We first estimated the physical conditions from the emission 
detected in the CO 3-2, 5-4, and 6-5 transitions both in HVC 
and LVC . We modelled each velocity component as a sim- 
ple uniform slab, adopting the size (filling factor) estimated 
above. Calculations were done in the large-velocity gradient 
approach, using the CO colUsional coefficients determined by 
Flower (2001) for ortho-H2 collisions in the range 5 K-400 K. 
For temperatures beyond 400 K, the collisional coefficients were 
extrapolated adopting a temperature dependence of Vr/400 K. 
Both components appear to have the same gas column den- 
sity N(CO) ^ 10" cm-2 (Fig. 4). We found the high-velocity 
gas component to be unambiguously associated with hot gas 
(T > 350 K) of moderate density x 3.0 x 10"* cm"^, whereas 
the low-velocity component arises from gas at a lower temper- 
ature (100 K) and higher density (~ 3.0 X 10^ cm"'). The tem- 
perature estimated for the extended component agree reasonably 
well with other determinations from NH3 and CH3CN (Tafalla 
& Bachiller 1995; CodeUa et al. 2009). 

With the physical conditions derived from the CO analy- 
sis, we modelled the integrated intensity and the fine profile 
of the 0-H2O lio - loi transition as well as the reported PACS- 
measured 179 fim H2O line intensity (~ 10"'^W cm"^, Nisini et 
al. 2010) to compute the total ortho water abundance in each 
velocity component. We used a radiative transfer code in the 
large- velocity gradient approach (and slab geometry) detailed in 
Melnick et al. (2008), taking into account an ortho to para ra- 
tio (OPR) of 1.2 for H2, as derived from Spitzer (Neufeld et al.. 



2009). Here, we assume the absorption component at -1-2.9 km 
is due to foreground gas unrelated to L1157-B1. Together, the 
two components of the 0-H2O 1 10 - loi fine produce a total 
H2O 2i2 - loi 179 Aim line flux of 1.1 x 10-'''Wcm-2. For 
the temperature range derived from our CO analysis, we esti- 
mated ortho-H20 column densities of (4.0 -5.0) x lO'^ cm"^ and 
(2.5 - 3.0) X 10'*^ cm-2 for the LVC and the HVC, respectively. 
Assuming an OPR of 3, we derived the H2O abundance from 
comparison with the gas column densities estimated from CO 
(see Table 1). We obtained an abundance ratio [H20]/[C0] ^ 
0.8 in the high-velocity gas, which is consistent with previous 
results from ODIN (Benedettini et al. 2002) and agrees reason- 
ably well with the predictions of steady-state C-shock models for 
this set of physical parameters (shock velocity ^ 20kms"', 
pre-shock density n(H2) = 5x10^ cm~^; Gusdorf et al. 2008a,b). 

An interesting prediction of our model is that the HVC con- 
tribution to the 179 yum flux dominates over the LVC contribution 
(see last column in Table 1). The higher temperature of this com- 
ponent drives the neutral-neutral reactions that efficiently form 
H2O, and the higher shock velocity can more efficiently remove 
water from grain mantles (see Melnick et al. 2008), resulting in 
the much greater ortho-H20 column density than in the LVC. 
Comparison with NH3 also suggests that the water production in 
the HVC is strongly dominated by high-temperature reactions 
(see CodeUa et al.). The higher ortho-H20 column density is 
what produces the higher 179 jxm line flux from this component. 
Consistent results are obtained by Nisini et al. (2010) based on 
a 179 yum PACS map and previous ODIN and SWAS observa- 
tions of the 0-H2O lio - Iqi line, assuming one single physical 
component dominates the water line emission in the HIFI beam. 
Follow-up observations of the higher-excitation lines of CO and 
H2O with HIFI will help us constrain more accurately the physi- 
cal conditions of each velocity component (density, temperature) 
and more generally in the shock. 
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